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The formula for the initial mass spectrum of primordial black holes (PBHs), which can be used 
for a general case of the scale dependent spectral index, and for a wide class of models of the 
gravitational collapse, is derived. The derivation is based on the Press and Schechter formalism. 
The comparative analysis of different types of initial mass spectra used in concrete calculations 
is carried out. It is shown that densities of background radiations (u, 7) from PBH evaporations 
depend rather strongly on a type of the gravitational collapse and on a taking into account the 
spread of horizon masses at which PBHs can form. Constraints on parameters of the primordial 
density perturbation amplitudes based on PBH evaporation processes and on atmospheric and solar 
■ neutrino data are obtained. 
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Studies of cosmological and astrophysical effects of primordial black holes (PBHs) are important because they enable 
one to constrain the spectrum of density fluctuations in the early Universe. If the PBHs form directly from primordial 
' density fluctuations then they provide a sensitive probe of the primordial power spectrum on small scales, > 10~ 9 pc. 

In particular, limits on PBHs production can be used to constrain models of inflation, in which the perturbation 
' amplitudes are relatively large at small and medium scales. 

In the simplest case, if we assume that the cosmological PBH formation is dominated, approximately, by primordial 
perturbations of one particular scale (i.e., there is some characteristic epoch of the PBH formation), we can obtain 
I limits on the initial mass fraction of PBHs, 

O' Pi = pPBH,i/ptot,i, (1-1) 

<n : 

where ppbh,i an d Ptot.% arc the PBH and total energy densities, respectively, at the time ti of the formation. This 
fraction can be expressed by the integral 

O ; ft = / P(6)dS, (1.2) 
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where P(S) is a probability distribution for density fluctuations entering horizon at ti, S is a density contrast, and 5 C 
is a minimum value of S required for the collapse. If the probability distribution is assumed to be Gaussian, one has 

P(S) = —^e~^ , (1.3) 
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where a is the rms fluctuation amplitude on a given scale. It is just this value that is determined by the primordial 
power spectrum. The connection between (3 and a is very simple in a case of the Gaussian distribution: 

ft«ae~i^. (1.4) 

The limits on ft arise from the entropy production constraints pi , from a distortion of the microwave background 
, from the cosmological nucleosynthesis constraints Q| . In these cases PBHs which give constraints have evaporated 
completely to the present time. On the contrary, the gravitational constraint (flpBH.o = Ppbh.o/Pc < 1), 



M 



1/2 



Pi, is valid for PBH masses M > I0 15 g which survive today. 

One should stress that all these limits are based just on the approximation that all PBHs form at the same scale 
and, correspondingly, the initial mass spectrum of PBHs is <5-function-like or, at least, is "non-extended" one. Only in 
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this case one can approximately express the observational constraints through the initial mass fraction Evidently, 
a more accurate treatment should operate with the initial PBH mass spectrum directly. 

The most strong limit on a PBH formation in the early Universe is due to the possible contribution of evaporating 
PBHs to the cxtragalactic 7-ray and neutrino backgrounds at energies ~ 100 MeV. The limit of such kind was 
obtained in the work of Page and Hawking p] , authors of which assumed that the differential initial mass spectrum 
of PBHs has power law form predicted in Carr's work ||, 

n B H(M BH ) = (a-2) f^f^) M^prQpBHM , (1.6) 
M* < M BH . 

Here, a — 2.5 for PBHs formed in radiation-dominated era, M* is the mass of a black hole whose life-time is equal to 
the present age of the Universe. 

This PBH mass spectrum is, clearly, the example of an "extended" spectrum (it was derived in || by considering 
the PBH formation as a process stretched in time). Naturally, the famous Page-Hawking constraint was formulated 
in terms of a PBH number density rather than in terms of Pi. According to la], the upper limit on the present PBH 
number density is ~ 10 4 pc~ 3 (or, that is the same, flpsH.o 10~ 8 ). This constraint was improved in later works [Q, 
where the same initial P BH mass spectrum was used as input. 



The derivation of Eq.(l.£) was based on the assumption of exact scale invariance (scale- independence of the per- 
turbation amplitudes). The Carr's work |jj appeared before an advent of the inflation hypothesis, and at that time 
it seemed improbable that the possible case of a growth of the perturbation amplitudes with a decrease of the scale 
can be of any importance (from the point of view of observational evidences of the PBHs existence) . Now we know 
that, due to inflation, minimum values of PBH mass in the initial mass spectrum can be rather large (10 13 — 10 14 g 
or even more). It means that the PBHs can have rather extended mass spectrum even if the primordial fluctuations 
are not strictly scale invariant. 

In general, initial PBH mass spectrum depends on cosmological and astrophysical aspects of the model used for 
its derivation (and, through the model, on such parameters as a time of the end of inflation, a reheating temper- 
ature, a spectral index of the density perturbations (or parameters of the inflationary potential) and, last but not 
least, on parameters, characterizing the process of the gravitational collapse leading to the PBH's birth). Evidently, 
observational constraints on the PBH production (especially those derived from measurements of extragalactic diffuse 
backgrounds of 7-rays and neutrino) can be used as constraints on at least one of these parameters. Clearly, these 
constraints depend on the initial PBH mass spectrum, i.e., on the parameters of this spectrum which are considered 
as free in a course of the constraint's derivation. 

Usually, the parameters characterizing primordial density fluctuation amplitudes (e.g., the spectral index, if it is 
scale independent) are objects of the constraining. Parameters of the gravitational collapse and a method used for a 
summation over epochs of the PBH formation (if such a summation is performed) determine, by definition, the type 
of the initial PBH mass spectrum. All other parameters (characterizing, in particular, the end of inflation and the 
beginning of radiation dominated era) are " external" and considered as free ones. 

In the present work we obtain some cosmological constraints following from the possible contribution of PBH evap- 
orations in extragalactic diffuse neutrino background. Assuming that the power spectrum of primordial fluctuations 
has a power-law form, P(k) ~ k n , we present the observational limits as constraints on values of the spectral index n 
(using the normalization on COBE data). In a more general case of the non-power P(£;)-dependence one can directly 
constrain parameters of a concrete inflation model. We consider, as an example, a case of the running mass inflation 
model Q and obtain constraints on parameters of the corresponding inflationary potential. 

The main physical assumption used in the work is that a PBH formation process occurs during radiation dominated 
era only. The main "external" parameter is ti, a time of the beginning of radiation era. At this moment of time we 
have, by assumption, the primordial spectrum of density perturbations, and no primordial black holes. In a short 
period between the end of inflation and ti there can be additional amplification of density perturbation amplitudes 
(e.g., in a preheating phase) but, by assumption, in this period there is no PBH formation. The time U is connected 
by a usual way with a corresponding initial temperature of radiation era, which we call a reheating temperature Trh- 

The paper is organized as follows. In Sec.|l] we derive, using the Press and Schechter formalism, the general formula 
for the initial PBH spectrum which is valid for any law of P(/c)-dependence and for a wide class of models of the 



gravitational collapse. In Sec. Ill we giv e the comparative analysis of different types of initial PBH mass spectra 



used in concrete calculations. In Sec. IV the neutrino background spectra from PBH evaporations are derived and 
the corresponding constraints on cosmological parameters are obtained. Discussions and conclusions are presented in 
Sec.0. 
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II. INITIAL PBH MASS SPECTRUM FORMULA 



According to the Press - Schechter theory 0, the mass distribution function n(M,5 c ), which is defined such that 
n(M, 5 c )dM is the comoving number density of gravitationally bound objects in the mass range (M,M + dM), is 
determined by the equations (see, e.g., ]10||) 



n(M,5 c ) = ^ 



dF 

8M^ 



(2.1) 



F(M,6 C ) = 



a R (M) 



(2.2) 



Here, <tr(M) is the rms mass fluctuation on a mass scale M, i.e., a standard deviation of the density contrast of the 
density field smoothed on a size R having mean mass M, F(M, 5 C ) is a fraction of the volume occupied by the regions 
that will eventually collapse into bound objects with masses larger than M (S c is a minimum value of the density 
contrast needed for developing a nonlinear growth of the density fluctuations in the overdense region). The density 
Pi is, in our case, the background energy density at the initial moment of time, t = ti (it is a moment of a beginning 
of the growth of the densit y flu ctuations) . 

By the construction, Eq^^l) supposes a step-by-step formation of the whole mass spectrum of gravitationaly bound 
objects (PBHs, in our case). Heavier PBHs form at later times when corresponding masses of overdense regions cross 
horizon. The nonzero probability that small PBHs are included in larger ones at later times is approximately taken 
into account by the diffe rentiation of F(M, 5 C ) with M. The correct normalization of the spectrum is ensured by the 
extra factor 2 in Eq.(2.2). 

Introducing the double differential distribution n(M, S) [pi] defined by the expression 



n(M, 8 C ) = I n(M, S)dS 



So 



one obtains from Eqs. fl2.1| - [2.2| ) 

n(M, 6) 



2 Pi 



nM <T r (M) 



da R (M) 



dM 



The rms mass fluctuation a R (M) is given by the integral p|] 

4W = / (Jgf 5Uk)W*(kR)T\k)^, 



(2.3) 



(2.4) 



(2.5) 



where W(kR) is the smoothing window function, Snik) is the (nonsmoothed) horizon crossing amplitude, T{k) is the 
transfer function taking into account the microphysical processing of density perturbations after entering horizon, a 
and H are a cosmic scale factor and an expansion rate. 

There are two time-dependent factors in the right hand side of Eq.(2.5): the factor {aH) A in the denominator 
describing a growth of the perturbations with time and the transfer function. We need the amplitude a R ( M) at the 
initial moment of time, t = ti. A role of the transfer function is essential at large k, if an integrand in Eq.( |2.5| ) is not 
cut off by a fc-dependence of the window function and <5y(fc). In the case when a top- hat window function is used 
fc™^ 1 with n > 1, the contribution of large k modes in the integrand is rather large and a convergence of 



and 5jj 



the integral in Eq.(2.5) is provided just by the transfer function. The form of this function at an epoch of the PBH 
formation is, of course, model dependent. One can assume, for simplicity, that at t = ti the transfer function has a 
simplest form: 



T(k)\ t=u =Q(k e -k), 



(2.6) 



where Q(x) is the step function, and the border value k e is of the order of k en d, the horizon scale a t the end of inflation 
|P^| . From our point of view, however, one cannot exclude the possibility that k e » k enc i- If Eq.( |2.6|) is assumed, the 
equation (|2.5|) can be rewritten, for a power-law of the primordial power spectrum and a top-hat window function, as 



4(M) 



M 
(aHY 



C 2 



n ■ 6 H( k fl) 



(2.7) 
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c 2 



, n 



x n+2 W 2 (x)dx = 9 / x n+2 



sm x cos x 



dx. 



(2.8) 



If the spectral index n is slightly more than one, n = 1.2, and if k e /k en d = 10 one has C ~ 5, i.e., the value of C is 
close to the corresponding value Q for the present epoch. Having no reliable model for the transfer function at early 
postinflationary epoch, we assume everywhere below, where a power-law form of P(k) is used, that C is constant, 
and equal to 4.7, throughout all times. Correspondingly, 



4(M) 



k JL- i CHl{k sl ) = 1 ^ I a 2 H {M h ) 



(aH) 



(any 



(2.9) 



In Eq^iO) k/i is the comoving wave number, characterizing the perturbed region, 



k f i = -, 



R 



(2.10) 



and Mh is the fluctuation mass at the moment when the perturbed region crosses horizon in radiation dominated era. 
The value of Mh is expressed through kfi (see below). The function an(Mh) introduced by Eq.(2.9) is, evidently, the 
smoothed horizon crossing amplitude (a standard deviation of the density contrast at a moment when the fluctuation 
crosses horizon). 

The general formula for af 1 (Mh), which is valid for an arbitrary ^-dependence of P(k), is 



x 3 W 2 (x)S 2 H {k f ix)dx. 



(2.11) 



At the moment th when the fluctuation crosses horizon one has 

Ra(t h ) w th, 

and, in radiation dominated era (a ~ yi), 



M h 



1 

K fl 



(2.12) 



(2.13) 



At the same time, an initial value of the fluctuation mass, M, is proportional to R 3 , so Mh ~ M 2 ! 3 . The exact 
connection is 



M h = Ml /3 M 2/3 

where Mi is the horizon mass at the initial moment of time, t = t*. 



(2.14) 



M t - U 



(diHi) 



2 ' 



The denominator (aH) in Eq.( |2 1 9j ) must be taken also at t = ti and, using the proportionality 



M ~ R 3 a 3 (U)p(U) ~ ^= ~ k n 3 -L, 



one can easily see that 



M 



{aiHi} 2 \M, 



-2/3 



Due to Eq. (|2.17|) , one has the simple connection between <jr and an- 



M 



-2/3 



(2.15) 



(2.16) 



(2.17) 



(2.18) 
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For the following we parametrize the fc-dependence of 5 B (k) by the formula 

n(k)-l 



Introducing now the new variable 5' , 



5 H (k) = 8 H (k ) — 
\k 



5' = 8 



M 
Mi 



2/3 



and a general form of the connection between PBH mass, fluctuation mass and density contrast, 

M BH = f(M,S'), 

one gets the PBH mass spectrum ti B h{Mbh)'- 



n B H{M BH ) = J n BH (M BH ,5')d5' , 



(2.19) 



(2.20) 



(2.21) 



(2.22) 



n BH {M BH ,5') =n(M,6) — - 



dS' dM bh 

Using Eqs.( 2.4), (2.9), (2.14), ( 2"7Ti| - 2.21), one obtains the final expression for the PBH mass spectrum 

dk fl 



nBHiMBH) = ^ pi lM^imr){^ 



—ln^- + n ( fc /') ~ 1 1 
2 U k 2 fy, 



Here the following notation is used: 



o\{M h ) 
! dn(k) 



<9M 



d/(M, <y')/dM 



The limits of integration in Eq.( [2.24 ) are determined from Eqs.( 2.2fj| ) and (2.21) using the conditions 



S m in = S c , M mm = M™ m = M l 



The connection between kfi and M is determined by the expressions 

fe/i = (a eq H eq ) [ -jj^T- 



aeqH eq — V2HoQ m Cl r ^ 2 , M e q = Q~T~t eq 



(2.23) 



(2.24) 



(2.25) 



(2.26) 



(2.27) 



together with Eq.(2.14). In these formulae a e(? and £/eg are the cosmic scale factor and the expansion rate at the 
moment of matter-radiation equality, t = t eq , Hq is the Hubble constant (Hq = lOOhkms^ 1 Mpc" 1 ) , M eq is the horizon 
mass at t eq , fi m and Q r are matter and radiation densities in units of p c . We consider in this paper, for simplicity, 
the case of flat models with zero cosmological constant. 



III. DIFFERENT TYPES OF INITIAL PBH MASS SPECTRA 



A concrete form of the initial mass spectrum of primordial black holes is determined, first of all, by peculiarities of 
the gravitational collapse near a threshold of the black hole formation. According to analytic calculations of seventies 
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| |T5| . pl| a critical size of the density contrast needed for the PBH formation, 5 C , is about 1/3. Besides, it was argued 
that all PBHs have mass roughly equal to the horizon mass at a moment of the formation, independently of the 
perturbation size. The approximate connection between the PBH mass and the horizon mass in such models is very 
simple @0: 

M BH = 1 l ' 2 M h , (3.1) 

where 7 is the ratio of the pressure to energy density (7 = 1/3 in radiation dominated era). 

However, it was shown recently that near a threshold of the black hole formation the gravitational collapse behaves 
as a critical phenomenon |jlq|. I n this case the initial mass function will be quite different from the analogous function 
in the analytic calculations [fL5| . The main feature is that the PBH mass formed depends on a size of the fluctuation 

B 

M BH = K M h (5' -6 c y\ (3.2) 

In this formula Mh is the horizon mass at the time when the fluctuation enters horizon; S c , n, 7 are parameters of a 
concrete model of the critical collapse p^ , p^| . The corresponding f(M, 8') function defined in Eq.(2.21) is determined 



using the relation fl2.14|) 



f(M, 8') = kMI /3 M 2 / 3 (8' - 5 c y*. (3.3) 

It is seen from Eq.( |3.2|) that the PBH mass may be arbitrarily small independently of the value of Mh- Besides, a 
value of the critical overdensity, 5 Cl in such models is typically ~ 0.7, i.e., about a factor of 2 larger than the value 
used in the standard collapse case. 

The second important ingredient of a PBH initial mass spectrum calculation is a taking into account the spread 
of horizon masses at which PBHs are formed. This problem exists independently of a nature of the gravitational 
collapse. 

As is stated in the Introduction, it is assumed very often that a majority of the PBH formation processes occurs 
at the shortest possible scale. In particular, the authors of ref. Jl(| determined the PBH initial mass function under 
the assumption that all PBHs form at the same horizon mass. The accuracy of such an approximation was studied 
in the work [ fl9| using the excursion set formalism, and it was found that it is good enough, at least in a case of 
the power-law density perturbation spectra with the spectral index close to 1.2 (i.e., slightly "blue" spectra, with the 
value of n satisfying observational constraints). 

As it follows from the previous Section, the effect of an accumulation of PBHs formed at all epochs after the 
beginning of ordinary radiation era can be conveniently taken into account by the Press-Schechter formalism . For 
a particular case of the critical collapse, the expression for the PBH mass spectrum, based on the Press-Schechter 
formalism, was obtained in refs. pp|,pi|. The spectrum formula is pi| 



n BH {M BH ) 



1/2^-5/2 



BH 



cr H (M h ) 



\'-2 



(3.4) 



An upper limit of the integration in Eq.(3.4) is determined by the expression 



(M BH V'\ 



1 



(3.5) 



The formula (3A) was derived using the assumption that the primordial power spect rum P{k) has a power-law 
form (i.e., the spectral index n is constant throughout all scales). Correspondingly, Eq.( |3.4| ) can be obtained from 
the general expression ( [2.24 ) omitting there the n'-terimThe amplitude a B (Mh) can be normalized on COBE data 
using the connection between a B and S B given by Eq.(|2.9|), 



a H (M h ) = CS H (k f i), 



(3.6) 



G 



and the relation P,pH 



5 H (k fl )^2-W- 5 



- 5 (M, 



< 'i 



Mh 
M, 



eg 



(3.7) 



Here Myo is the present horizon mass. 

All numerical calculations the results of which are presented in this paper are carried out with the following values 
of the parameters [M : 



Ik = 0.36 



0.7. 



(3.8) 



The known formula |U| for the PBH mass spectrum in the standard collapse case can be obtained from Eq.(^J) 
by using the substitutions plfl 



Ik 







7 



1/2 



7 



and the approximate relation 



S' 2 

ds' | V - 1 



s' 2 



Substituting Eqs. (|3.9|) and ( p.lOl) in Eq.(|3.4|) one obtains 

+ 3 [2 



nBH(M BH ) 



^ Pi Ml/ 2 M- 5 J 2 a^(M h ) 



exp 



2a%(M h 



(3.9) 



(3.10) 



(3.H) 



The Carr's formula [Eq.(|l.6|)] follows from Eq.( 3.1l| ) if two additional assumptions are used: i) the amplitude 
an(Mh) does not depend on Mh (exact scale invariance) and ii) there is no minimum value of PBH mass in the initial 

5/2 

spectrum. In this case the spectrum is AM BI j , as it follows from Eq.(3.11). The factor A is found by normalization 
of the present PBH mass spectrum on the present value of ppbh,o = Qpbh,oPc, using the approximate formula for 
the instantaneous PBH spectrum at t = to ll]], 



n BH (m,t c 



(Af3 + m 3 ) 2 / 3 
(O(x) is the step function), and the expression 



A{Ml + to 3 )- 5/2 9(to - A/*), 



(3.12) 



oo 

Ppbhm = J mnBH(rn,t )dm. 



(3.13) 



We carried out in this wo rk t h e calc ulations of the neutrino background spectra from PBH evaporations using the 
mass spectra given by Eqs.(3.4),( 3.11 ) and, in parallel, the PBH mass spectrum from ref. fill , 



n BH (Mbh) = 



p .(MMM> lh " 
y kMi 



\{-!TM—) +*<= 



2TTa H (M h )M B HMijk 



(3.14) 



In this expression the amplitude <th is determined by the same formula as above, but with the constant value 

of Mh, Mh — Mi, in accordance with the assumption of au t hors of ref. [[16] that all PBHs form at the smallest scale. 

We have, for calculations of ubh(Mbh), using Eqs.(3.4),( 3.11 ) and ( 3.14 ), two free parameters: the spectral index 
n, giving the perturbation amplitude through the normalization on COBE data on large scales, and ti, the moment 
of time just after reheating, from which the process of PBH formation started. The value of ti is connected, in our 
approach, with a value of the reheating temperature, 



0.301 9 : 



-1/2 M p i 
T 2 

RH 



(3.15) 
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~ 100 is a number of degrees of freedom in the early Universe). The initial horizon mass Mj is expressed through ti 



Mi Si -Trtfpi. (3.16) 



Typical results of PBH mass spectrum calculations, for the case of a scale independent spectral index, are shown on 
Fig.l. We use, for convenience, the following abbreviatio ns: K L, NJ and BK signify, corres p ondin gly, t he PBH mass 
spectra calculated by formulae of Kim and Lee |Q (Eq.( 3.11 )), Niemeyer and Jedamzik Q (Eq.( 3.14 )) and Bugaev 



and Konishchev |TT|] (Eq.(^J)). One can see from Fig.l that at large values of the spectral index (n = 1.3) there is a 
clear difference in a behavior of two mass spectra, NJ and BK, based both on a picture of the critical collapse: BK 
spectrum is less steep at Mbh > Mi. Of course, this difference is not so spectacular in cases when the spectral index 
n is more close to 1. 

The common feature of NJ and BK spectra is the maximum at Mbh — Mi. Besides, both spectra have the same 
slope in the region Mbh < Mi. The distinctive peculiarity of K L spectra is a sharp cut off at Mbh — Mi, due to the 



rigid connection between PBH mass and horizon mass [Eq.( 3.1)1 m the standard collapse case 



IV. NEUTRINO BACKGROUND SPECTRA FROM PBHS EVAPORATIONS 

Evolution of a PBH mass spectrum due to the evaporation leads to the approximate expression for this spectrum 
at any moment of time: 

2 

ns "( m >*) = (3cd + m3)V3 nBH ( i3at + m3 ) 1/3 ) ' ( 4J ) 

where a accounts for the degrees of freedom of evaporated particles and, strictly speaking, is a function of a running 
value of the PBH mass m. In all our numerical calculations we use the approximation 

a = const = a(M B l H x ), (4.2) 

where Mg H x is the value of Mbh in the initial mass spectrum corresponding to a maximum of this spectrum. Special 
study shows that errors connected with such approximation are rather small. 
The expression for a spectrum of the background radiation is |Tl[| 

3 r ™2 



- (i) / "" v»:t,3)3/3 "°" [<m + - 3 ) 1 ' 3 ] ■ w ■ (i + .),»>.-**■» 

= J F(E,z)dlog w (z + l). 



(4.3) 



In this formula dj, a and do are cosmic scale factors at ti, t and at present time, respectively, and f(E,m) is a 
total instantaneous spectrum of the background radiation (neutrinos or photons) from the black hole evaporation. It 
includes the pure Hawking term and contributions from fragmentations of evaporated quar ks an d from decays of pions 
and muons (see jll| and earlier papers f22|-p4|| for details). The exponential factor in Eq.([D|) takes into account an 
absorption of the radiation during its propagation in space. The processes of the neutrino absorption are considered, 
in a given context, in re f. pif . 

In the last line of Eq. ( |4.3| ) we changed the variable t on z using the flat model with J7a = for which 

(4.4) 

fi r = (2.4 • itfh 2 )- 1 , h = 0.67. 
On first five figures (Figs.|[0) some results of calculations for the case of a scale independent spectral index are 



presented. Several examples of z-distributions (the integrands of the integral over z in Eq.(4.3)) are shown on Fig. 2. 
Again one can see the rather strong difference of two cases: the existence of a tail of heavy masses in BK spectrum 
leads to a relative enhancement of low z-contributions. The effect becomes more distinct with a rise of the reheating 
temperature. 
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The sharp cut-off of all ^-distributions near z ~ 10 7 is entirely due to the neutrino absorption The shrinkage of 
z-distributions at large Trh in NJ-case is due to the absence of large masses in the spectrum (PBHs of small masses 
evaporated earlier, and their radiation today is more redshifted). 

Some results of calculations of the neutrino background spectra are shown on Figs. 3-5. The functional form of these 
spectra is qualitatively the same as in the standard collapse case pl|: E~ 3 law at large neutrino energies and a flat 
part at low energies with the crossover energy depending on Trh- The only new feature is the sharp steepening of 
the spectra at high values of Trh in NJ case (Fig. 4) connected with the corresponding shrinkage of the z-range. 

Fig. 5 shows the sensitivity of background neutrino intensities to a chosen value of the spectral index n, at a hxed 
Trh- One can see from this figure that the neutrino background flux from PBH evaporations is comparable with the 
theoretical curve |2j| for the atmospheric neutrino flux at ~ 100 MeV if n ~ 1.29 (at Trh ~ 10 9 GeV). Such large 
values of the spectral index are, probably, excluded by the recent large-scale experiments. 

Next figure shows the constraints on the spectral index following from data of neutrino experiments, as a function 
of the reheating temperature. As in our previous work, we used for a description of these constraints the data of the 
Kamiokande atmospheric neutrino experiment [Eq] and the experiment on a search of an antineutrino flux from the 



Sun [|27| (see (ll) for details) 



For completeness, we show on the same figure the constraints following from extragalactic diffuse gamma ray data 



part of them (NJ case) is in qualitative agreement with the results of ref. 28 1. It is seen that the constraints are much 
more weak in NJ case (at Trh > 10 10 GeV). It is again the result of an absence of the tail of large masses in NJ 
spectrum, leading to a shrinkage of the z-distributions and to a steepening of the background neutrino spectra. 

One can see from Fig that, in general, spectral index constraints following from the comparison of neutrino 
background predictions with existing data of neutrino experiments are rather weak if the purely power law of primordial 
density fluctuations is assumed for all scales (and, correspondingly, the normalization on COBE data is used). However, 
it does not mean that neutrino background intensities from PBH's evaporations cannot be noticeable. If, for example, 
the spectrum of primordial density fluctuations is a combination of two simple power law spectra pQ], i.e., 



<r H (M h ) ~ M h 4 , M h < M hc < M eq ; 



(4.5) 



a H (M h )~M h ' 1 , M eq >M h >M hc , 

then the small value of n; ~ 1 (ni -lwO) which follows, in particular, from COBE data, does not contradict with 
the possibility that n s — 1 is large and, correspondingly, an at small scales is also large. If Mh c and n; are known 
from some model, one can easily obtain from the neutrino data the constraints on the small scale spectral index 
ri s , in the same manner as it is done above for the all scale spectral index n. Several cosmological models of such 
kind, which predict a large density fluctuation amplitude just on small scales (while the amplitude on large scales is 
constrained by the small CMBR anisotropy) appeared recently. The general problem of all these models is that the 
large density perturbations on small scales can lead to an overproduction of the primordial black holes. An origin of 
the large density perturbations and the scale where they are close to maximum depend on the model. For example, 
hybrid inflation models with two stages of inflation p9[ ] predict a characteristic spike in the perturbation spectrum, 
and the model parameters (those determining a duration of the second stage of inflation) can be chosen so that PBHs 
evaporations give large contributions to extragalactic backgrounds of neutrinos and 7-quanta. Inflation models with 
a subsequent preheating phase also predict perturbation spectra growing strongly toward small scales ]30|] (as a result 
of the resonant enhancement of field fluctuations during preheating). 

At the end of this section we consider, with more details, the case of running mass inflation models, which also give 
a spectral index with rather strong scale dependence. The inflationary potential in these models is of the form |Sy: 



V = V 



, 1 < 

1 _ 77 c " 
2 I 



/,2 



(4.6) 



where c and Q are positive constants, <p is inflation field. This form corresponds to a loop correction coming from 
softly broken global supersymmetry. The constant c is, essentially, the coupling strength of the field in the loop which 
is expected to be of the order of 0.01 — 0.1 |pl| . The spectral index is given by the simple expression: 

~ [n{k) - 1] = ae- cN W - c , (4.7) 

where N(k) is a number of e - folds of inflation that occur after the epoch when a scale k leaves horizon outside. The 



new parameter a in Eq.(4.7) is expressed through Q, c and N(k). Using the connection d\nk = —dN one obtains 

!„' = ™ e -c*m. (4.8) 
2 k x ' 
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For normalization of the perturbation amplitude to COBE data we use |3!2J the parameter Ncobe = 
lrt(k enc [/kcoBE) connected with N(k) by the expression 

N(k) = Ncobe - In . k , (4.9) 

KCOBE 

where the COBE scale (the center of the range explored by COBE) is 

kcoBE = ^hMpcr 1 . (4.10) 
We can fix the parameter Nc obe a ssuming the instant reheating after the end of slow-roll inflation, and expressing 



Ncobe through T RH using Eq.(2.27) 



k ?r = (a eq H eq )^L, (4.11) 



Mi ' 

the connection between Mi and Trh [Eqs.( [3T5| ),( |3.16| )l, and the relation 

kjr = k end = e N °° BE a H . (4.12) 

Instead of this, we put, for simplicity, Ncobe = 45 (for all possible values of Trh), i-e., we assume that, in a general 
case, kj\ ax < k enc i (radiation era begins not immediately after the end of slow-roll inflation). The slow-roll condition, 

V = lM^<l, (4.13) 

holds if values of the parameter c are not too large. If, e.g., ucobe = 1, this condition can be rewritten as 

7] w ce cNcOBE < 1. (4.14) 



The horizon crossing amplitude 8n{k) is calculated using Eq.(2.19), where we put fco = kcoBE- The expression (4.7) 
for n{k) has two parameters: a and c. If we suppose that the spectral index at k = kcoBE is known from experiment 
(e.g., according to ref. p3| , it follows from COBE data that ucobe — 1.02 ± 0.24), we exclude one parameter (e.g., 
a). Now we can use the same steps as above to obtain the constraints on the p aram eter c. The only differe nce i s 
that now one must use the general expressions for the PBH mass spectrum, Eq.( [2.24 ), and for &jj(Mh), Eq.( [2.11 ) 



These constraints are shown on Fig.^ for three possible values of tlcobe- Solid lines on the figure correspond to the 
critical collapse case, dotted line (calculated for the value ucobe — 1-0 only) correspond to the standard picture of 
the collapse. The lines show the upper limits of the parameter c. One can clearly see the difference between two 
cases: the sharp rise of the dotted curve at Trh ~ 10 8 GeV is directly connected with the sharp cut-off of the initial 
mass spectrum (for the standard collapse case) at Mbh = Mi. Indeed, due to the proportionality 

M i ~ *i ~ T RH, ( 4 - 15 ) 

at small enough reheating temperatures the minimum values of the PBH mass in the spectrum are too large and, 
therefore, the large part of PBHs cannot evaporate until the present time. 

Finally, on Fig. |] we show the resulting constraint curves for the horizon crossing amplitude Sn^k), for fixed values 
of ucobe and Trh (for the critical collapse case). The curves correspond to the upper limits of the par amet er c 



shown on the previous figure. The end points of the power spectrum on the figure is determined from Eq.(4.11). It 
is seen from Fig.|| that allowed values of 5^ at the shortest scale (for a given Trh) are rather large, 

S 2 H {k?i ax ) ~3-10" 4 . (4.16) 

The corresponding a H value is - 0.08. At T RH = 10 10 GeV the shortest scale, k™ ax - 10 18 {h Mpc" 1 ) - k endl 
corresponds to a rather small value of the horizon mass, M™ n = Mi ~ 10 11 g. PBHs of such a small mass evaporated 
at t - 10 6 s. 
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V. CONCLUSIONS 



The main conclusion of the work is the following: constraints on parameters of cosmological models from evapo- 
rations of primordial black holes depend rather strongly on a form of their initial mass spectrum. Number densities 
of PBHs predicted by the models are extremely sensitive to values of primordial density perturbation amplitudes 
and, therefore, even very steeply falling mass spectra of PBHs can give useful constraints. As it follows from FigJ|, 
constraints on the scale independent spectral index n at high values of a reheating temperature are drastically dif- 
ferent in NJ and BK cases. We thoroughly traced the origin of this result: the absence of a long high mass tail 
in the initial NJ spectrum (Fig.|l|) leads to a relative shrinkage of the redshift distributions of evaporated neutrinos 
(Fig.^J) and, as a consequence, to a steepening of the neutrino background spectra and to a decrease of the background 
intensities (Fig.|j), resulting, eventually, in a weakening of the spectral index constraint. So, the constraint values at 
Trh > 10 10 GeV on Fig.O are entirely due to an effect of the summation over all epochs of the PBH formation in the 



approach of ref. |20 21j] . Analogously, we see, on an example of the parameter c of running mass inflation models, 
that the constraints at low values of Trh also depend on a type of the gravitational collapse: in the standard collapse 
case the constraint at Trh < 10 8 is practically absent (dotted curve on Fig.fj]) due to a corresponding absence of a 
low mass tail in the initial PBH mass spectrum of this case (compare KL and BK spectra on Fig.[l]; the behavior of 
the corresponding PBH mass spectra in a case of a scale dependent n is qualitatively similar). 

The initial number density of PBHs is larger in a case of the standard picture of the gravitational collapse (as 
compared with the critical collapse case). It is due to relatively small value of the critical overdensity (5f = 7 = 1/3). 
In the region near the maximum of PBH mass spectra the ratio of intensities is given by the approximate relation 

(KL) 3g- T 2 

"/;// 



(BK) 



e 2 "h . (5.1) 



n BH 



Correspondingly, intensities of the neutrino and gamma backgrounds produced by PBH's evaporations are smaller 
and spectral index constraints are systematically weaker in the critical collapse case (see, for the comparison with the 
standard collapse case, Fig. 10 of ref. jT^j). 

It is worth noting that spectral index constraints followed from neutrino evaporations are stronger at high Trh 
values (as compared with the constraints from 7-quanta), especially in the NJ case (Fig.ra. The reason is simple: 
at high Trh the relative contribution of large redshifts in z-distributions is very large (FigH) and, correspondingly, 
the absorption of v, 7 during propagation in space is important. Naturally, absorption effects are more sufficient 
for 7-quanta than for neutrinos (for 7-quanta z max ~ 700, while for neutrinos z max ~ 10 7 JO], as we can see, in 
particular, from curves of Fig.||). 

One should note that, in general, the constraints on a scale independent spectral index followed from PBH evapo- 
rations (Fig||) are rather weak. Probably, so large deviations of the spectral index from 1 (n ~ 1.28) are excluded by 
the latest data (see, e.g., pl|). Therefore, constraints from PBH evaporations are more useful for those cosmological 
models in which the spectral index is scale dependent. In such models large density perturbation amplitude at small 
scales (and, correspondingly, large effects from evaporations) can coexist with the small value of Sn(k) at COBE scale 
(Fig.||). We showed in this paper, using, as an example, the running mass inflation model, that in this case the PBH 
evaporation process can be used for a constraining of model parameters [in particular, parameters of an inflationary 
potential, (Fig^)]. We stress, once more, that all constraints of such type depend on assumptions used in a derivation 
of the initial PBH mass spectrum. 
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FIG. 1. Examples of PBH mass spectra for two values of Trh- Solid lines: KL spectra [18], dotted lines: NJ spectra [16], 
dashed lines: BK spectra [11]; n = 1.30 in all cases. 
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FIG. 2. Red-shift distributions (integrands of the expression (4.3) for the background spectrum) for neutrino energy 
E — 1CP 1 GeV (solid lines) and E = 10~ 2 GeV (dashed lines), for PBH mass spectrum (left column) and BK mass spec- 
trum (right column), for three values of Trh', n = 1.30 in all cases. The dimension of F(z) is (s~ 1 cm~ 2 sr~ 1 GeV~ 1 ). 



14 



- 20 

> 

O 




-4 -3-2-10 1 2 

log 10 (E/1 GeV) 



FIG. 3. Neutrino background spectra for three values of Trh, for BK mass spectrum; n = 1.30. 
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FIG. 5. Neutrino background spectra for three values of the parameter n, for BK mass spectrum. Dashed line represents the 
theoretical atmospheric neutrino spectrum for Kamiokande site (averaged over all directions) [25] . 
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FIG. 6. Constraints on the spectral index n as a function of the reheating temperature Trh- Solid and dotted lines: BK 
and NJ mass spectra, atmospheric u e and solar v e experiments. Dashed and dot-dashed lines: BK and NJ mass spectra, 
extragalactic diffuse gamma-ray background data. 
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FIG. 7. Upper limits on values of the parameter c of the running mass inflation model, for different reheating temperatures 
and for several possible values of tlcobe- Solid lines: critical collapse case, dotted line: standard collapse case. 
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FIG. 8. Constraints on the power spectrum of primordial density fluctuations in the running mass inflation model, for several 
values of ticobe and for T RH = 10 10 GeV. All curves are calculated for the critical collapse case. 
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